*replication of the tables in "Labor Market Information and Parental Attitudes toward Women Working Outside the Home: Experimental Evidence from Rural Pakistan" as of August 3, 2022

*table1
*using wwoh_pak.dta
sort No
by No: gen hid=_n
local vars f_edu m_edu f_age m_age nhh kammee farmland v_livestock v_homeplusland fac_corr1 fac_lc flfpRF 
sum `vars' if Year2==0&hid==1&Treated==1
sum `vars' if Year2==0&hid==1&Treated==0
gen Tg=Treated*Tgroup
gen Ti=0
replace Ti=1 if Treated==1&v_treatment_type==2
gen clst=villageNo if Tg==1|Vpercent==0|Vpercent==1
replace clst=HHID+40 if clst==.
***balance test following Heß (2017)& 3rd comment***
ritest Treated _b[Treated], reps(1000) cluster(clst): reg f_edu Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg m_edu Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg f_age Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg m_age Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg nhh Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg kammee Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg farmland Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg v_livestock Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg v_homeplusland Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg fac_corr1 Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg fac_lc Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg flfpRF Treated if Year2==0&hid==1, cluster(clst) robust

**using roster_baseline.dta
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==0
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==1
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==0
***balance test following Heß (2017)& 3rd comment***
gen Tg=Treated*Tgroup
gen clst=villageNo if Tg==1|Vpercent==0|Vpercent==1
replace clst=HHID+40 if clst==.
ritest Treated _b[Treated], reps(1000) cluster(clst): reg hh_age_ind Treated if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg hh_enroll_ind Treated if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg hh_edu_ind Treated if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&hh_enroll_ind!=1, cluster(clst) robust

*using wwoh_pak.dta
gen glfp3=0
replace glfp3=1 if occ==7|occ==8
gen temp0=glfp3 if Year2==0
gen temp1=glfp3 if Year2==1
sort HHID hhindid
by HHID hhindid: egen glfp3_0=total(temp0)
by HHID hhindid: egen glfp3_1=total(temp1)
drop temp0 temp1
gen glfp4=0
replace glfp4=1 if occ==3|occ==7|occ==8
gen temp0=glfp4 if Year2==0
gen temp1=glfp4 if Year2==1
sort HHID hhindid
by HHID hhindid: egen glfp4_0=total(temp0)
by HHID hhindid: egen glfp4_1=total(temp1)
drop temp0 temp1
*table1 outcome variables
local vars info5 info6 info7 flfppre 
sum `vars' if Year2==0&hid==1&Treated==1
sum `vars' if Year2==0&hid==1&Treated==0
sum glfp3 glfp4 if Year2==0&Treated==1&hh_age_ind>=15&hh_age_ind<30
sum glfp3 glfp4 if Year2==0&Treated==0&hh_age_ind>=15&hh_age_ind<30
***balance test following Heß (2017)& 3rd comment***
ritest Treated _b[Treated], reps(1000) cluster(clst): reg info5 Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg info6 Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg info7 Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg flfppre Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg glfp3 Treated if Year2==0&hid==1, cluster(clst) robust
ritest Treated _b[Treated], reps(1000) cluster(clst): reg glfp4 Treated if Year2==0&hid==1, cluster(clst) robust

*table2 Panel A following 2nd comment (and 3rd comment)
reg info5_re1 Treated if Year2==0&hid==1, cluster(clst)
reg info6_re1 Treated if Year2==0&hid==1, cluster(clst)
reg info7_re1 Treated if Year2==0&hid==1, cluster(clst)
reg flfppre_1 Treated if Year2==0&hid==1, cluster(clst)
reg glfp3_1 Treated if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)
*multihypothesis testing
rwolf2 (reg info5_re1 Treated if Year2==0&hid==1, cluster(clst)) (reg info6_re1 Treated if Year2==0&hid==1, cluster(clst)) (reg info7_re1 Treated if Year2==0&hid==1, cluster(clst)) (reg flfppre_1 Treated if Year2==0&hid==1, cluster(clst)) (reg glfp3_1 Treated if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)), indepvars(Treated, Treated, Treated, Treated, Treated) usevalid cluster(clst) reps(1000)

*table2 Panel B following 2nd comment
reg info5_re1 Treated info5_re0 if Year2==0&hid==1, cluster(clst)
reg info6_re1 Treated info6_re0 if Year2==0&hid==1, cluster(clst)
reg info7_re1 Treated info7_re0 if Year2==0&hid==1, cluster(clst)
reg flfppre_1 Treated flfppre_0 if Year2==0&hid==1, cluster(clst)
reg glfp3_1 Treated glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)
rwolf2 (reg info5_re1 Treated info5_re0 if Year2==0&hid==1, cluster(clst)) (reg info6_re1 Treated info6_re0 if Year2==0&hid==1, cluster(clst)) (reg info7_re1 Treated info7_re0 if Year2==0&hid==1, cluster(clst)) (reg flfppre_1 Treated flfppre_0 if Year2==0&hid==1, cluster(clst)) (reg glfp3_1 Treated glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)), indepvars(Treated, Treated, Treated, Treated, Treated) usevalid cluster(clst) reps(1000)

*table2 Panel C following 2nd comment
gen f_edu2=0 if f_edu==1|f_edu==2
replace f_edu2=1 if f_edu==3
replace f_edu2=2 if f_edu==4|f_edu==5
replace f_edu2=3 if f_edu==6|f_edu==7
gen m_edu2=0 if m_edu==1|m_edu==2
replace m_edu2=1 if m_edu==3
replace m_edu2=2 if m_edu==4|m_edu==5|m_edu==6
local vars i.f_edu2 i.m_edu2 f_age m_age nhh kammee farmland ls residence fac_corr1 fac_lc flfpRF
reg info5_re1 Treated info5_re0 `vars' if Year2==0&hid==1, cluster(clst)
reg info6_re1 Treated info6_re0 `vars' if Year2==0&hid==1, cluster(clst)
reg info7_re1 Treated info7_re0 `vars' if Year2==0&hid==1, cluster(clst)
reg flfppre_1 Treated flfppre_0 `vars' if Year2==0&hid==1, cluster(clst)
reg glfp3_1 Treated glfp3_0 hh_age_ind `vars' if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)
rwolf2 (reg info5_re1 Treated info5_re0 `vars' if Year2==0&hid==1, cluster(clst)) (reg info6_re1 Treated info6_re0 `vars' if Year2==0&hid==1, cluster(clst)) (reg info7_re1 Treated info7_re0 `vars' if Year2==0&hid==1, cluster(clst)) (reg flfppre_1 Treated flfppre_0 `vars' if Year2==0&hid==1, cluster(clst)) (reg glfp3_1 Treated glfp3_0 hh_age_ind `vars' if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(clst)), indepvars(Treated, Treated, Treated, Treated, Treated) usevalid cluster(clst) reps(1000)

*table2 Panel D following 2nd comment
local vars intT Year2 Treated i.f_edu2 i.m_edu2 f_age m_age nhh kammee farmland ls residence fac_corr1 fac_lc flfpRF
reg info5 `vars' if hid==1, cluster(clst)
reg info6 `vars' if hid==1, cluster(clst)
reg info7 `vars' if hid==1, cluster(clst)
reg flfppre `vars' if hid==1, cluster(clst)
reg glfp3 `vars' hh_age_ind if hh_age_ind>=15&hh_age_ind<30, cluster(clst)
rwolf2 (reg info5 `vars' if hid==1, cluster(clst)) (reg info6 `vars' if hid==1, cluster(clst)) (reg info7 `vars' if hid==1, cluster(clst)) (reg flfppre `vars' if hid==1, cluster(clst)) (reg glfp3 `vars' hh_age_ind if hh_age_ind>=15&hh_age_ind<30, cluster(clst)), indepvars(intT, intT, intT, intT, intT) usevalid cluster(clst) reps(1000)

*table3
reg info5_re1 intTsex Treated Tsex info5_re0 if Year2==0&hid==1, cluster(villageNo)
lincom intTsex+Treated
reg info6_re1 intTsex Treated Tsex info6_re0 if Year2==0&hid==1, cluster(villageNo)
lincom intTsex+Treated
reg info7_re1 intTsex Treated Tsex info7_re0 if Year2==0&hid==1, cluster(villageNo)
lincom intTsex+Treated
reg flfppre_1 intTsex Treated Tsex flfppre_0 if Year2==0&hid==1, cluster(villageNo)
lincom intTsex+Treated
reg glfp3_1 intTsex Treated Tsex glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(villageNo)
lincom intTsex+Treated
rwolf2 (reg info5_re1 intTsex Treated Tsex info5_re0 if Year2==0&hid==1, cluster(villageNo)) (reg info6_re1 intTsex Treated Tsex info6_re0 if Year2==0&hid==1, cluster(villageNo)) (reg info7_re1 intTsex Treated Tsex info7_re0 if Year2==0&hid==1, cluster(villageNo)) (reg flfppre_1 intTsex Treated Tsex flfppre_0 if Year2==0&hid==1, cluster(villageNo)) (reg glfp3_1 intTsex Treated Tsex glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(villageNo)), indepvars(intTsex, intTsex, intTsex, intTsex) usevalid cluster(villageNo) reps(1000)
rwolf2 (reg info5_re1 intTsex Treated Tsex info5_re0 if Year2==0&hid==1, cluster(villageNo)) (reg info6_re1 intTsex Treated Tsex info6_re0 if Year2==0&hid==1, cluster(villageNo)) (reg info7_re1 intTsex Treated Tsex info7_re0 if Year2==0&hid==1, cluster(villageNo)) (reg flfppre_1 intTsex Treated Tsex flfppre_0 if Year2==0&hid==1, cluster(villageNo)) (reg glfp3_1 intTsex Treated Tsex glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(villageNo)), indepvars(Treated, Treated, Treated, Treated) usevalid cluster(villageNo) reps(1000)

*table4
reg info5_re1 Tg Ti info5_re0 if Year2==0&hid==1, cluster(villageNo)
test Tg=Ti
reg info6_re1 Tg Ti info6_re0 if Year2==0&hid==1, cluster(villageNo)
test Tg=Ti
reg info7_re1 Tg Ti info7_re0 if Year2==0&hid==1, cluster(villageNo)
test Tg=Ti
reg flfppre_1 Tg Ti flfppre_0 if Year2==0&hid==1, cluster(villageNo)
test Tg=Ti
reg glfp3_1 Tg Ti glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(villageNo)
test Tg=Ti

*tableA1
sum bicycle motorbike car washmachine sewmachine generator tv ac cellphone fridge

*tableA2
sum o3 o4 o5
tab o3
tab o4
tab o5

*tableA3
local vars s_flfp1 s_flfp2 s_flfp3 s_flfp4
sum `vars' if Year2==0&hid==1&Treated==1
sum `vars' if Year2==0&hid==1&Treated==0
ritest Treated _b[Treated], reps(1000) strata(villageNo): reg s_flfp1 Treated if Year2==0&hid==1, a(villageNo) robust
ritest Treated _b[Treated], reps(1000) strata(villageNo): reg s_flfp2 Treated if Year2==0&hid==1, a(villageNo) robust
ritest Treated _b[Treated], reps(1000) strata(villageNo): reg s_flfp3 Treated if Year2==0&hid==1, a(villageNo) robust
ritest Treated _b[Treated], reps(1000) strata(villageNo): reg s_flfp4 Treated if Year2==0&hid==1, a(villageNo) robust

***Appendix R&R3
***balance test father, mother treated for table3***
*tableA4
local vars f_edu m_edu f_age m_age nhh kammee farmland v_livestock v_homeplusland fac_corr1 fac_lc flfpRF 
sum `vars' if Year2==0&hid==1&Treated==1&Tsex==1
sum `vars' if Year2==0&hid==1&Treated==1&Tsex==0
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg f_edu Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg m_edu Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg f_age Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg m_age Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg nhh Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg kammee Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg farmland Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg v_livestock Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg v_homeplusland Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg fac_corr1 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg fac_lc Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg flfpRF Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
***outcome variables balance test father, mother treated***
local vars info5 info6 info7 flfppre 
sum `vars' if Year2==0&hid==1&Treated==1&Tsex==1
sum `vars' if Year2==0&hid==1&Treated==1&Tsex==0
sum glfp3 glfp4 if Year2==0&Treated==1&hh_age_ind>=15&hh_age_ind<30&Tsex==1
sum glfp3 glfp4 if Year2==0&Treated==1&hh_age_ind>=15&hh_age_ind<30&Tsex==0
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg info5 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg info6 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg info7 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg flfppre Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg glfp3 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg glfp4 Tsex if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
***using roster_baseline.dta balance test father, mother treated***
*tableA4
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tsex==1
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tsex==0
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tsex==1
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tsex==0
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg hh_age_ind Tsex if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg hh_enroll_ind Tsex if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1, cluster(villageNo) robust
ritest Tsex _b[Tsex], reps(1000) cluster(villageNo): reg hh_edu_ind Tsex if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&hh_enroll_ind!=1&Treated==1, cluster(villageNo) robust

***balance test group, individual treated for table 4
*tableA5
local vars f_edu m_edu f_age m_age nhh kammee farmland v_livestock v_homeplusland fac_corr1 fac_lc flfpRF 
sum `vars' if Year2==0&hid==1&Treated==1&Tgroup==1
sum `vars' if Year2==0&hid==1&Treated==1&Tgroup==0
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg f_edu Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg m_edu Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg f_age Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg m_age Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg nhh Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg kammee Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg farmland Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg v_livestock Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg v_homeplusland Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg fac_corr1 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg fac_lc Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg flfpRF Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
***outcome variables balance test group, individual treated***
local vars info5 info6 info7 flfppre 
sum `vars' if Year2==0&hid==1&Treated==1&Tgroup==1
sum `vars' if Year2==0&hid==1&Treated==1&Tgroup==0
sum glfp3 glfp4 if Year2==0&Treated==1&hh_age_ind>=15&hh_age_ind<30&Tgroup==1
sum glfp3 glfp4 if Year2==0&Treated==1&hh_age_ind>=15&hh_age_ind<30&Tgroup==0
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg info5 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg info6 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg info7 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg flfppre Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg glfp3 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg glfp4 Tgroup if Year2==0&hid==1&Treated==1, cluster(villageNo) robust
***using roster_baseline.dta balance test group, individual treated***
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tgroup==1
sum hh_age_ind hh_enroll_ind if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tgroup==0
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tgroup==1
sum hh_edu_ind if hh_relate_ind==4&hh_enroll_ind!=1&hh_age_ind>=15&hh_age_ind<30&Treated==1&Tgroup==0
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg hh_age_ind Tgroup if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg hh_enroll_ind Tgroup if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&Treated==1, cluster(villageNo) robust
ritest Tgroup _b[Tgroup], reps(1000) cluster(villageNo): reg hh_edu_ind Tgroup if hh_relate_ind==4&hh_age_ind>=15&hh_age_ind<30&hh_enroll_ind!=1&Treated==1, cluster(villageNo) robust

*tableA6
gen T1=0 
replace T1=1 if Vpercent==.4
gen T2=0
replace T2=1 if Vpercent==.6
gen T3=0
replace T3=1 if Vpercent==1
gen Treated1=T1*Treated
gen Treated2=T2*Treated
local vars T1 T2 T3 Treated1 Treated2
reg info5_re1 `vars' info5_re0 if Year2==0&hid==1, cluster(villageNo)
reg info6_re1 `vars' info6_re0 if Year2==0&hid==1, cluster(villageNo)
reg info7_re1 `vars' info7_re0 if Year2==0&hid==1, cluster(villageNo)
reg flfppre_1 `vars' flfppre_0 if Year2==0&hid==1, cluster(villageNo) 
reg glfp3_1 `vars' glfp3_0 if Year2==0&hh_age_ind>=15&hh_age_ind<30, cluster(villageNo)
